?7a(1) Anomaly at high temperature: the scalar-pseudoscalar splitting in QCD 
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We estimate the splitting between the spatial correlation lengths in the scalar and pseudoscalar 
channels in QCD at high temperature. The splitting is due to the contribution of the instanton/anti- 
instanton chains in the thermal ensemble, even though instanton contributions to thermodynamic 
quantities are suppressed. The splitting vanishes at asymptotically high temperatures as AM/M oc 
(Aqcd /T) b , where b is the beta function coefficient. 
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There has been a lot of interest lately in the high temperature phase of Quantum Chromodynamics [1]. Although 
the main thrust of the recent activity, triggered by the RHIC data, has been at temperatures close to the transition 
. temperatures, there still are open questions about the high temperature deconfined phase of QCD. This regime 
may soon be experimentally accessible; the recently published RHIC data [2] suggests a temperature in the range 
O ,. 300 < T < 600 Me V, considerably above the expected transition temperature of T c w 170 MeV. At asymptotically 
high temperatures the basic QCD physics is perturbative, however some aspects of it cannot be understood without 
nonperturbative effects. The generation of magnetic mass is an example of such an effect [3]. Another, simpler, set 
of questions has to do with contribution of instantons at high temperature [4] . Although the instanton contribution 
i— i . to thermodynamic quantities like pressure and energy is negligible, they must give the leading contribution to certain 
correlation functions. In particular, the question of whether and when the axial Ua{1) symmetry of QCD is restored 
above the phase transition has been debated for some time. Clearly, since the anomaly is an operator equation, it 
O i' remains true on the operator level also at nonzero temperature [5]. At low temperatures, the relation between the 
anomaly equation and ir° decay has been elucidated in [6]. At very high temperatures, the effects of the anomaly 
are however generally believed to be very small because of increasing irrelevance of instanton contributions at high 
temperature [4]. A convenient measure of the explicit breaking of the Ua(1) symmetry in the theory with two massless 
flavors is the difference of the spatial correlation lengths in the scalar and pseudoscalar channels. This quantity has 
been calculated on the lattice and has been found to be small and probably nonvanishing [7-10], although the lattice 
studies are rather inconclusive because the signal is numerically delicate, especially in the presence of chiral zero 
modes. We are, moreover, unaware of an analytic calculation of this splitting. The aim of this short note is to provide 
t-H ' such a calculation and to demonstrate that the instanton configurations contribute decisively to this quantity even 
, though their contribution to pressure etc. may be negligible. This work is inspired by [11] where a very similar 
calculation was performed in the three dimensional Georgi-Glashow model, which also exhibits a deconfinement phase 
transition [12]. Indeed there is a very close analogy between the symmetry structure of the 3D Georgi-Glashow model 
and QCD with massless quarks, as discussed in detail in [13]. 
, We are interested in calculating the equilibrium (equal time) correlation functions 
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(S(x)S(y)) , (P(x)P(y)) (1) 

where 

S(x)=$ i 1> i = $il>? + ji?tf ; P{x)=$ i iwl> i = il$:tf-<$?tf) (2) 
The flavor index i takes values 1, 2. The action of the axial UaO-) symmetry as usual is 

tf^etotf ; Vf^e--Vf (3) 
If the axial anomaly is absent, the diagonal correlators must vanish in the axially symmetric thermal ensemble 

(#(a#f '(x) tf(y)^(y)) = o = {$?{x)tf{x) $%)#(y)> (4) 

and therefore 

(S(x)S(y)) = (P(x)P(y)) (5) 

Although the axial symmetry is anomalous, the effects of the anomaly are "transmitted" to the QCD spectrum 
via the nonperturbative instanton and anti-instanton contributions [14-17]. At zero temperature instantons are finite 
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action solutions of Euclidean equations of motion. In the regular gauge the vector potential of an instanton in the 
SU(2) gauge theory is 

W = 5^ (6) 

where the totally antisymmetric 't Hooft symbol is given by r] a ^ = e afJiI ^; p, v = 1,2,3; rj afl 4 = S afJi . The instanton 
in the SU(N) theory is essentially the same configuration embedded into the SU(N) group. All color orientations of 
the instanton have to be summed over in the path integral. 

Due to conformal symmetry of the classical QCD action, the instanton solution is characterized by a size parameter 
p. The action of an instanton is independent of p classically 

S,-£ (7, 

The most important quantum correction to this result stems from the fact that the coupling constant runs with scale, 
and so semiclassically 

s - - m < 8 » 

The measure of integration over the instanton size is dn oc ^(pA) b , where b = (11N C — 2Nf)/3 is the coefficient of 
the one loop /3-function. The measure is peaked towards large size instantons, which makes instanton calculations at 
zero temperature uncontrollable. At zero temperature therefore, although we understand qualitatively that the mass 
splitting between the scalars and pseudoscalars is due to instanton-like fluctuations in the vacuum, a first-principles 
controllable calculation is not possible. The best one can do at present is the semi phenomenological instanton liquid 
model [18, 19], which works quite well. 

At finite temperature the situation is similar, as long as the temperature does not significantly exceed the non- 
perturbative QCD scale Aqcd [17]. However, at high temperature the instanton calculations are much more under 
control since the temperature provides an effective infrared cutoff on the instanton size [20] . In the imaginary time 
formalism at finite temperature instantons are periodic solutions of the classical equations of motion with period 
(3 = l/T in the imaginary time direction. Instantons with core size p <C j3 look essentially the same as at zero tem- 
perature, but large instantons with p » (3 look very different. The core of size (3 is accompanied by a dyon-like field 
in the spatial dimensions, and so at large distances a large size instanton looks like a monopole [4]. Most importantly, 
since the nonzero temperature T provides an external scale, already classically the instanton action depends on the 
instanton size. As a result, large instantons with sizes p > (3 are exponentially suppressed. Since the contribution to 
the path integral due to very small instantons is suppressed by the zero temperature measure, this means that the 
main contribution to physical observables comes from the instantons with the size p ~ f3. 

An instanton at zero and finite temperature is accompanied by fermionic zero modes. There is one left-handed 
zero mode for ipi of each flavor and one right-handed zero mode for ipi of each flavor. Due to the presence of these 
zero modes every tunneling event associated with the instanton is accompanied by the change of the axial JZa(1) 
charge in the vacuum, giving a concrete manifestation of the anomaly. The presence of the zero modes also modifies 
significantly the interaction between the instantons and anti-instantons. In pure Yang-Mills theory, an instanton (I) 
and an anti-instanton (A) interact weakly with the interaction "potential" decreasing as a power of the distance. 
Put differently, a field configuration corresponding to an Instanton/ Anti-instanton pair separated by a distance R 
(assuming the sizes of the instantons are much smaller than the separation) has the action 



S IA (R) ~ ™£ 

g 2 



2 2 

1 + f lA 



R A 



(9) 



The actual formula is somewhat more complicated since the instantons also posess a color orientation, and thus the 
interaction depends on the relative color orientation of the instanton and an anti-instanton. We disregard this subtlety 
since at high temperature the favored orientation is such that I and A are parallel in the color space [17, 19]. 

In the presence of the fermion zero modes the I-A interaction becomes long range. At zero temperature the wave 
function of the zero mode decreases as 1/R 3 away from the center of the instanton, and as a result the I-A interaction 
is logarithmic. At high temperatures on the other hand, the fermionic zero mode wave function has a characteristic 
exponential decay in space with the '"mass"' equal to the lowest Matsubara frequency ttT. As a result the I-A 
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interaction potential becomes linear. For an I-A pair separated by a large distance R in the spatial direction the 
interaction potential at high temperature due to 27V/ zero modes is 

S IA (R)~ S I + S A + 2N f irTR (10) 

At these high temperatures the instantons are like monopoles, and the "confining" linear potential between a monopole 
and an anti-monopole binds the monopoles into pairs of net zero topological charge. It is believed that the transition 
to the "molecular" phase coincides with the chiral symmetry restoration and the deconfinement phase transition and 
happens around T c « 170 Mev in the theory with two massless flavors [17]. 

A simpler model which exhibits a very similar behavior and instanton binding above the deconfinement phase 
transition is the Georgi-Glashow model in 2+1 dimensions [11-13]. Although in this model the instantons are bound 
in pairs at high temperature, they nevertheless have a direct effect on some observables. As was shown in [11], in 
the 2+1 dimensional Georgi-Glashow model they give the main contribution to the mass (spatial correlation length) 
splitting between the scalar and the pseudoscalar channels. We expect similar effects in QCD. Instanton effects should 
be there in a theory with any number of flavors, but they are simplest to see in the two flavor case, so in the following 
we therefore concentrate on Nf = 2. 

Consider the equal time correlation function 

g L l(x - y) = ®?{x)tf(x) ^(y)^(y)) (11) 

As mentioned above, in the absense of the Ua{^) anomaly this correlation function vanishes, since it has a net axial 
Ua(1) charge. Thus perturbative calculation of this correlator at high temeperature gives a vanishing result. It is 
however easy to see that this correlator does not vanish on a single Instanton configuration. The instanton has four 
zero modes - exactly the right number to be saturated by the fermionic operators in the observables. Thus, as long 
the points x and y are very far away from each other, \x — y\ <C 1/T, an instanton centered at a point a gives a 
contribution 

<#(*)Vf (*) (y))i « aT e e- s 'e- 2 ^ S - S ^y-^ (12) 

where a is a constant of order one. The maximal contribution comes from the instantons located along the straight 
line connecting points x and y. The instanton contribution does not depend on the position of the instanton along 
the line. Integrating over the position of the instanton with the weight Tdz we obtain 

$?(x)tf(x) i>*{y)tf{y))i « (aT 6 )(cT)\x-y\e- s 'We- 2 " T ^ (13) 

where we have set the instanton size to be equal to the inverse temperature and have allowed for a dimensionless 
constant c, of order one, which must appear due to the integration over the instanton sizes. It should be possible to 
estimate c using the dependence of the instanton action on p [4] , but our main interest in this note is to establish the 
parametric dependence of the correlation length on temperature. The next contribution comes from a configuration 
of two instantons and one anti-instanton, alternating along the straight line connecting x and y. The zero modes 
are again completely saturated and so this contribution does not vanish. The main contribution comes from the 
configuration where two of the zero modes of each instanton are saturated by the zero modes of the anti-instanton. 
The typical distances between the instantons clearly is greater than their size, and thus we can neglect all interactions 
between them except for that induced by the fermionic zero modes. The contribution of this configuration gives 

$?(v)tf{y))iAi « {aT^-c^x-yfe-^e- 2 ^-^ (14) 
The leading contribution from a configuration with an alternating chain of n instantons and n — 1 anti-instantons is 

$?{x)tf{x) $?{y)ti{y)) InA n-i « (aT^cTe- 8 ^) 2 "-^- 2 ^-^ f dy n f" dx n . x [" ,U Sl , ... H d Vl 

J X J X J X J X 

(15) 

Summing all the contributions we obtain 

(4>*(x)tf(x) 4>?{y)tf(y)) « (aT e )e~ 2 * T ^ sinh [cTe~ s ^\x - y\) (16) 

As already mentioned, we have assumed that all the instantons and anti-instantons have the same orientation in color 
space, because at high T the integration over the orientations is dominated by the configuration where instantons and 
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anti-instantons are parallel [17]. The integration over orientations will only affect the value of the pre-exponential 
factor c. 

Turning our attention to the axially symmetric correlator, we note that perturbatively 

(tf (x)^(x) ^f(y)^(2/)>perturbative « («T 6 ) (17) 

The origin of the exponential factor is now not the zero modes but rather just the contribution of the free lowest 
Matsubara frequency contribution. The constant factor aT 6 is the same as in eq.(12) etc. since it originates from 
the disconnected piece associated with points x and y separately. Clearly instanton/anti-instanton chains give a 
contribution to this correlator in exactly the same way as before. The only difference is that now the number of 
instantons and anti-instantons in the chain must be equal in order to saturate all the fermionic zero modes. The reult 
is then 

{tf(x)i/>?(x) i>?{y)tf{y)) « (aT e )e- 2 « T ^ cosh (cTe~ s '\x - y\) 
From equations (16) and (18) we find 

(S{x)S(y)) oc e - M sl*"fl ; (P(x)P(y)) K e -M*\*-v\ ; (S(x)P(y))=0 

with 

M s = 2ttT - cTe~ SlW ; M P = 2irT + cTe~ Sl{l3) 
Thus we find that the splitting between the scalar and pseudoscalar correlation lengths is given by 

M \ T J 

which is the advertised result. In the end, this is essentially the expected semiclassical quantum mechanical splitting 
determined by exponentiation of the single-instanton action [21, 22]. The key observation is that at high T the 
problem reduces to one of a linear chain of alternating instantons and anti-instantons, with the total number being 
cither different by one (anti-symmetric correlator), or equal (symmetric correlator). This is very similar to the behavior 
found in the Georgi-Glashow model [11], as we expected. 

We note that although we have discussed isosinglet correlators, the results apply also to isovectors since the chiral 
SU(2) (g> SU(2) symmetry is restored in the high temperature phase. Thus, for example we also have 

M aa - ( K Q cp \ b 

Although this ratio is not parametrically suppressed at T ~ 2T C , the large power b ~ 10 may explain the fact that at 
these temperatures the difference in the correlation lengths is difficult to detect in lattice calculations [10]. It would 
be very interesting if the high T regime could be probed further with high statistics lattice computations. 

A simple consistency check on our calculation is to add a 9 term to the QCD Lagrangian. In the chiral limit 
the spectrum of masses and correlation lengths should not depend on the value of 9 as it can be eliminated by the 
anomalous £7(1) rotation. This means that with 9 ^ we should recover the same correlation lengths, but in the 
correlators of the axially rotated operators. We can include the effect of the QCD 9 angle in our calculation by noting 
that the correlators with one extra instanton or anti-instanton in the chain will acquire a phase, while the correlators 
with equal numbers of instantons and anti-instantons will be unchanged. Thus, Gll acquires an extra phase e* e , Grr 
acquires an extra phase e , while Glr and Grl are unchanged. Then we find 

(S(x)S(y)) oc 2e- Mlx ~ vl [cosh(AM|x - y\) - cos (9 sinh(AM|a; - y\)] 
(P(x)P(y)) oc 2e- Mlx - yl [cosh(AM|x - y\) + cos (9 sinh(AM|a; - y\)\ 

(P(x)S(y)) = (S(x)P(y)) oc 2e~ Mlx ~ yl sin (9 sinh(AM|x - y\) (23) 

with M = \{M$ + Mp), AM = \{Mp — Ms). The eigenvalues of the correlator matrix are the same as before while 
the eigenvectors are rotated precisely by the axial rotation, as expected. 

We conclude with two brief comments. First, for a theory with more than two massless flavors, the I-A chains 
do not contribute to the correlation function of two fermionic bilincars, since the number of zero modes on an 
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instanton does not match the number of fermionic operators. Instead they contribute to the Green's functions with 
Nf bilinears. Second, at high temperature in QCD the instantons can be thought of as Skyrmions of the field A% 
in the dimensionally reduced theory [24]. This is very similar to the situation in the 3D Georgi-Glashow model, 
where instanton- monopolcs became vortices of the abelian part of A . The interesting difference is that it is easy 
to understand the linear interaction between vortices [11, 23] as a direct result of the nonexistence of a continuous 
U(l) symmetry. On the other hand it is more difficult to understand a linear potential between Skyrmions, since the 
Skyrmion field decays fast at infinity. Since this linear interaction is the consequence of the fermionic zero modes, it 
would be interesting to understand how this interaction arises in the effective Lagranian of the dimensionally reduced 
theory due to the integration of the quark fields. 
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